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Abstract 



A Bose-Einstein condensate may be used to make precise measurements of 
weak forces, utilizing the macroscopic occupation of a single quantum state. 
We present a scheme which uses a condensate in a double well potential to 
do this. The required initial state of the condensate is discussed, and the 
limitations on the sensitivity due to atom collisions and external coupling are 
analyzed. 
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I. INTRODUCTION 



The experimental demonstration of Bose-Einstein condensation (BEC) in trapped quan- 
tum gases of alkali atoms [|I]-|4]] opens up great opportunities in atomic physics, condensed 
matter physics and quantum optics. As a macroscopic quantum object, the physical proper- 
ties of the Bose-Einstein condensate (BEC) have naturally attracted great attention. While 
much experimental and theoretical work has already been done to investigate the basic 
properties of the condensates, possible practical application is still an open question. Very 
recently, experimentalists at JILA have attempted to build an atomic clock, based on per- 
sistent tunneling oscillations between two coupled BECs Although this type of atomic 
clock is still crude, the experiments have stimulated the search for practical applications of 
an atomic BEC. 

In this paper, we propose a scheme to detect weak forces by employing a BEC confined in 
a double- well potential. The basic idea for this purpose is shown in fig.(|l|). The condensate 
is prepared initially in a coherent superposition of the extremal eigenstates of the operator 
of particle number difference between the two wells. The requirement for such a state will 
be discussed below, but it should be noted that this state is an entangled state and cannot 
be described by semiclassical mean field theory. Under the action of a weak-force for a time 
r, the condensate will experience a phase shift. The phase shift can be detected by using the 
technique analogous to the Ramsey interference. The interference fringes can be read out by 
performing a homodyne measurement of the optical phase shift due to the dispersive inter- 
action of the condensate localized in one well and an optical field mode. We have analyzed 
the limitations on the accuracy of the scheme and show that a high-precision measurement 
can be achieved if the condensate contains a large number of coherently condensed atoms. 
This result is due to the quantum entanglement || inherent in the initial superposition state 
of BEC. 



II. TWO MODE SYSTEM 



We consider the case where a condensate has formed in a quartic double-well potential: 



V(x) = b (x 2 



(1) 



which has minima at x — ±x and a trap frequency of u = -y/86/mxo- In a two mode 
approximation where the total atom number N = N is conserved, the system hamiltonian 



can be described in terms of angular momentum operators ||,|10 

H = hQJ z + 2Tik.J 2 x . 
The commutation relations for these operators are 



J it Jj 



(2) 



(3) 



where allows cyclic permutations of i, j, k e {x, y, z}. 

The operator J x gives the condensate particle number difference between the two wells, 
J y corresponds to the momentum induced by tunneling, and J z is the difference in occupation 



2 



between the upper and lower energy eigenstates of the potential. The splitting between these 
levels is given by Q, which is the tunneling frequency, and k corresponds to the strength of 
the interparticle hard-sphere interactions. 

The two mode approximation is valid in the regime where the overlap between the single 
particle ground state modes of each well is small and where the many-body effects do not 
affect significantly the properties of these modes. These conditions lead to PJITfl 



— «1, N«£, (4) 
too H 

where a is the scattering length which determines the strength of the two-atom collisions 



and where r = J%/2mujQ characterizes the size of the wells. 



III. MEASUREMENT SCHEME 

As is shown in fig. ([I]), the detection of the weak force proceeds in several stages. The 
first step is to prepare the initial state of the system in a quantum state which can optimize 
the precision of the measurement by controlling the atomic parameters. With this in mind, 
we consider the weak tunneling limit Q << kN. In this case, the ground state of the 
hamiltonian (eq.(fj)) for attractive interactions (k < 0) is a superposition of the extremal 



eigenstates of the operator J x | 1 13] : 



\E) = -L(\j,-j) x + \j,j) x ). (5) 

We will see that with such an initial state the precision of the measurement can be controlled 
by the total number of atoms in the condensate. This state is easily represented on the Bloch 
sphere (fig. (|l](a))) by two diametrically opposed points on the equator. In principle, such 
a state could be prepared by allowing the atoms to condense into the ground state of a 
double-well potential. By increasing the height or width of the barrier, the tunneling rate 



could be adiabatically decreased [|13|] until the ground state evolves into the superposition 
state (eq.(|5|)). 

The second step in the measurement process, after having achieved the ground state 
described by eq.(^), and before the tunneling is turned on, is to allow the weak force to act 
for a certain time r. This force may be due to a varying gravitational field, for example, 
and has the effect of adding a linear ramp to the potential: 

H F = hAJ x , (6) 

where A is the frequency shift induced by the weak force (see fig.(|TJ(b))) and j = N/2. To 
avoid the self-phase shift due to interatomic collisions (nonlinear term in the hamiltonian 
eq.(|])), a technique such as that using Feshbach resonances |14| may be used to tune the 
size of the interatomic interaction close to zero at this stage. As a result, the initial ground 
state (eq.flBD) evolves into the superposition under the weak force: 

Mr)) = ± (e iA3> \j, -J) x + \j,j) x ) (7) 
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The state now contains the phase-shift induced by the weak force. This phase shift is 
"amplified" 2j times the single-particle value due to the specially prepared initial state. 

In the next step we propose to use a technique analogous to Ramsey interferometry. In 
order to detect the induced phase shift, we rotate the state (eq.(|7|)) around the Bloch sphere 
by 90°. Such a rotation can be achieved by turning on the tunneling between the two wells of 
the trap for a time t n /2 = If the interatomic interaction is strong, the rotation operation 
will be affected by the nonlinear term in eq.(|2[), which will distort the final state from the 
y-axis on the Bloch sphere and affect the precision of the measurement. To avoid this, we 
continue to use the Feshbach resonance to suppress the collisions. The state is then rotated 
to 



iAjr 



V2 



(8) 



Finally, a number measurement can be performed on one of the condensates. This 
corresponds to a projection onto the J x eigenstates. The resultant probability distribution 
is 



PJm) 



JArj 



(j,m\j, -j) v + e 



-iArj 



(9) 



Now the inner product of the J x and J y eigenstates will be a binomial function of m peaked 
around m — 0, with 



x(j,m\j, -j) 



T x(j,m\j,j) t 



(10) 



This leads to 



P x (m) = 2cos 2 (Ajr + mvr/2) x (j,m\j,j) 



2 cos 2 Ajr 
2 sin 2 Ajr 



A.i- ./•./) 

x(j,m\j } j) i 



m even 



m odd 



(11) 
(12) 



which describes how the output fringes are shifted by the presence of the force A. In the 
absence of the force, all the odd fringes are absent, but for A ^ 0, the probability that a 
particular measurement will fall on an odd fringe is 



Pr(odd) = sin 2 Ajr ~ (Ajr) 



(13) 



for small A. 



IV. MEASUREMENT READOUT 



The measurement of atom number is effected through a homodyne scheme ||10|| . The 
condensate is placed in an optical cavity, which at the time of the readout stage of the 
measurement contains a light field which is highly driven and damped. The optical field 
is thus in a coherent state with amplitude cuo- The light field is detuned from any atomic 
resonance and so the condensate merely imposes a phase shift on the light. This can be 
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detected by measuring the quadrature components of the field. For a dispersive interaction 
which acts on a timescale tj over which the dynamics of the condensate itself is negligible, 

Hi = xJx^a, (14) 

where X is the measurement strength and a', a are the light field creation and annihilation 
operators. The quadrature components X = + a and Y = i(a) — a) then execute simple 
harmonic motion: 

X(tj) = cos( X tiJ x )X(0) + sm( X tiJ x )Y(0) (15) 
Y(tj) = cos( X tiJ x )Y(0) - sm( X tiJ x )X(0). (16) 

After time tj, the light field is rapidly damped out. If all the light is emptied from 
the cavity, then the integrated photocurrent for the measured quadrature (X) is the true 
distribution for the quadrature [[H||16|]. In terms of the Wigner function W, this is given by 
the marginal distribution. 

I f°° 

p(x) = - dyW{x,y), (17) 

4 J—oo 

where x = a* + a and y = i(a* — a). 

In general, the Wigner function will be a sum of gaussians, weighted by the atom number 
distribution of the condensate: 

^(a,«*) = - f P x (m)e- 2|tt -< |2 (18) 

II m=-j 

where a' m = aoe~ lxtim = ^(x(0)+iy(0))e~ lxtim . For convenience, we set the initial conditions 
of the light field such that x(0) = 0. This gives 

9 

W(x,y) = - V p x ( m ) e -^ x - y ^ sinxt ' m ^-^ y - y ^ cosxtim ^. (19) 

n m=-j 

After integrating, the marginal distribution is 

p(x) = -^= P x (m)e-^ x - y{0)sinxtim)2 . (20) 

Thus each m value is mapped onto a gaussian at position y(0) sin^t/m with width equal to 
one. This is illustrated in fig.(|2]). To be able to distinguish without ambiguity different m 
values in the output, there should be at least 4 standard deviations between the means of 
adjacent gaussians. The resultant condition on the atom- light coupling is then 

4 2 

Xti > = ^ (21) 

\y( Q )\ Kl 

and 



Xh « ^- (22) 
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If the fringes close to m = only are needed, then this last condition may be relaxed 
considerably to a condition which merely prevents aliasing: 

Xti < ^- (23) 

Nevertheless, there is a limit on the size of the induced phase shift. 

The analysis above assumes perfect detector efficiencies and an infinite time of integration 
so that all the light is removed from the cavity. The results can be generalized to hold when 



this is not the case |]17| . For a detector efficiency of r/oo and a total integration time of T, 



then the distribution for the integrated photocurrent is 

p(x) = -= i , »e" 5(l "' w(0)si " xl " n)2/ '', (24) 

V 27Tr l m =-j 

where r\ = r/^ (l - e"^ T ) and 7 is the damping rate of the cavity. The lower bound on the 
atom-light coupling becomes 

4 2 

Xti > n . 7-r> = —, | . (25) 

Thus the effects of detector inefficiencies and a low damping rate can be overcome by starting 
with a large coherent state amplitude in the cavity. For example, with detector efficiency of 
?7oo — -5, N p = 10 7 photons in the cavity and a measurement strength of x — 10 _2 s , then 
for 7T ^> 1, the lower limit on the interaction time is 

ti > 90ms (26) 

This value of \ is calculated using a trap frequency of ujq/2-k = 32Hz, beam waist w = 
30/im, light detuning 5/2ir = 100MHz, saturation intensity I s = 17W/m 2 , optical frequency 
a;/27r = 3.8 x 10 l4 Hz, atomic linewidth r a /27r = 10 7 Hz and incident power P = 6mW, in 
a cavity 10cm long. 



V. PHASE ERRORS AND PHASE DIFFUSION 

Interatomic collisions are necessary to produce the initial superposition of condensates in 
the first stage of the scheme, but their effect on subsequent stages of the measurement scheme 
is unwanted. We have assumed that they can be suppressed using a Feshbach resonance. 
However, it may be unfeasible to use this technique, so we now discuss briefly the effect of 
the nonlinearity on the different stages of the scheme. 

In the second part of the measurement scheme, when the tunneling is turned off, the size 
of the self-interaction term may be comparable to that of the weak force term. This will 
induce an extra self-phase change: 

= -MA. (27) 
However, since this self-phase change is the same for both components of the superposition, 



there is no net effect on the output probability distribution (eq. fljjD ). 
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In the third stage, due to the constraints of the two-mode approximations, there are 
limits to the size of the interwell coupling. If the atomic collisions are weak compared to the 
tunneling term, then the nonlinear term will cause a collapse (through dephasing) in any 
tunneling oscillations. However, this should not occur before there is time for at least one 
quarter of an oscillation (a tc/2 pulse) to occur. If the self- interactions are stronger, then 
they will induce a nonlinear rotation around the J x axis and a diffusion of the distribution 
on the Bloch sphere. The effect of the extra rotation may be negated by adjusting the time 
of the pulse so that the final state lies in the J y — J z plane. The effect of the diffusion can 
not be so removed, and may wash out the interference fringes. 

Finally, there is last stage of the measurement when the condensate interacts with the 
cavity field. If there is no tunneling, the collisions will have no direct effect on the state of 
the light field, since the k J 2 term in the hamiltonian doesn't effect the x-distribution. If the 
Q is not exactly zero then, for a strong atom-light interaction, a back action may develop 
over time which will induce, through momentum fluctuations, tunneling ||10|| . This would 
directly affect the phase of the cavity field, and it also may open a way for the atom-atom 
interaction to have an effect. 

The effect of the back action on the condensate may be seen in the master equation for 
the system in which the dynamics of the optical field has been adiabatically eliminated [ITS 



p = -iQ[J z ,p) -i2K,[J x ,p] 
2y 2 \a\ 2 

+ i X \a\ 2 [J x1 p] - A 1 1 [J x ,[J x ,p]]. (28) 
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The last term in this equation is the decoherence induced by the external coupling (to the 
cavity field). Note that since it only involves J x operators, if Q is zero then the decoherence 
cannot affect the x-distribution and hence the induced light shift in the field. 

To avoid the waiting time involved in switching on the cavity field during the measure- 
ment process, it may be necessary to have the cavity on before the measurement begins. This 
will mean that the decoherence is active during the detection stage of the scheme, and so a 
random phase will be imparted to the condensate superposition. The double commutator 
in the master equation may be simulated by a stochastic term in the hamiltonian: 

b yff dt ^ V ' 

where dW is the Weiner increment. The resultant phase change has a standard deviation 
that grows with the square root of time: 

*^ = ^/T = T^a (0(r)) • (30) 

Thus the relative error caused by this phase diffusion may be minimized by increasing 
the detection time r or decreasing the strength of the interaction with the optical field. For 
typical parameters (as used above), with j = 50 and r = 160ms, the phase error is a$ = 0.08 
rad, which could be a major restriction on the sensitivity of the measurement. Hence it may 
be better to switch on the cavity only when it is time to use the optical field. 
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VI. MEAN FIELD LIMIT 



The scheme outlined above depended on starting in a state which was a quantum su- 
perposition of two condensates and on the resulting entanglement. For comparison, we now 
present the mean field analogue to show what features of this scheme remain in the absence 
of quantum entanglement. 

In the mean field limit, the system may be described by a Gross-Pitaevskii Equation 
MM (GPE): 

ih$(x,t)= (^-^^ + Rx + V(x) + U mx,t)\ 2 ^$(x,t), (31) 

where M is the atomic mass, the constant R is the gradient of the single particle potential 
due to the force and U is the strength of the interparticle interactions. Suppose that the 
interatomic collisions are negligible. Then, when the overlap between the wells is small, we 
may expand the mean field in terms of the local wavefunctions of each well: 

$(x,t) = &i(t)ui(a;) + b 2 {t)u 2 (x) (32) 

where 

Uj (x) = e -*M^ e (,-(-iF,o) 2 K ro = JIKL (33) 

n ' (2vrr )i V2Ma;o V ; 

and 

bj{t) = f u*Ax)$(x,t)dx. (34) 



The ground state energy of each of the local modes is E . From the GPE (eq.(|3T| ), the 
resultant equations of motion for the bj(t)s are 

6 i (t) = " ( " 1 ^ JtC ° 6 i (t) + yW*) (35) 

(36) 

As before, to perform the weak force measurement, we allow the force to act for a time 
r and then the tunneling for a time t n / 2 = 

bj ( T + U /2 ) = e-^^HjiO) cos ^ + i e -^ i+HRa: ^%^(0) sin ^ (37) 

(38) 



Suppose we start off with an equal occupation in each well, such that &i(0) = 62(0) = y N/2. 

Then the mean population difference is shifted by the presence of the weak force: 

(m) =^(\b 2 (r + t 7T/2 )\ 2 -\b 1 (r + t 7T/2 )\ 2 ) (39) 
N . 2Rx t 

= -- B n— (40) 
N 

= -—smAr. (41) 
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This did not occur in the previous quantum treatment, in which (m) = 0. Even when 
the system is simply in a number state (not a superposition) with an equal number of atoms 
in each well, i.e. \j, 0) x , the mean population difference is unaffected by the presence of 
the force. This is a demonstration of the fact that, as a classical treatment, the mean field 
situation cannot be regarded as the large number limit of a quantum number state. In 
quantum optics, it is the minimum uncertainty coherent state |a) which is most like a mean 
field with amplitude a. The analogue in this case is the atomic coherent state, or Bloch 
state: 



E 



2 j \ pm+J 

m + j J (l + \/3\ 2 y 



\j, m) 



2 ' 



(42) 



where j3 can be described in terms of the angular coordinates of a point on the sphere 
f3 = tmOe**. 

Consider the state given by /3 — 0, which is symmetric with respect to the two wells. This 
state is also the ground state \j, —j) z of the system when Q ^> \k\N, in other words, a state 
in which the coherence between the two condensates is well established through tunneling. 
If the system begins in this state, then after the measurement procedure, the difference in 
occupation between the two wells is as given above in the mean field approach (eq. ((4l|)). 

It may seem better to use this Bloch state as the initial state, since it may be easier to 
generate than the superposition state previously used (Eq.(|5|) and the sign of the weak force 
may be determined from the measurement of (m). However, the size of the induced phase 
change given in Eq.(^T|)is not amplified by j. In other words, the macroscopic occupation of 
a single condensed state in not being fully utilized. 

A comparison of the relative uncertainty in either case clearly demonstrates this point. 
In the first case, where the superposition state is used, the phase is inferred by the propor- 
tion of detection events falling on odd fringes (Eq.([l3|)). The uncertainty in this binomial 
distribution with probability P = sin 2 Ajr is 



5P 



! P(1-P) 
iVp ' 



(43) 



where Np is the number of detection events. The relative uncertainty in the phase is thus 

-l 



dP 



5P 



(44) 



When the initial state is the Bloch coherent state \ j, — j), the uncertainty in the mean of the 
distribution is 



S (m) = 

This gives a relative uncertainty in <fi of 

5(f) Id (m) 



\ 



(m 2 ) — {mY 

iVp 



(45) 



5 (m) 



1 



(46) 
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Thus the precision of the measurement grows in proportion to the number of condensed 
atoms when the entangled state (Eq.@) is used, but only as the square root of the number 
of atoms when the coherent state \ j, —j) is used. This demonstrates the advantage of using 
quantum entanglement of two macroscopically distinct states in order to make a highly 
sensitive force detector. 



VII. CONCLUSIONS 

If the weak force in question is gravity, then, with the parameters quoted above, the size 
of the induced phase shift is 

- = — - ~ 2.4 x lOVad per g, 47 

9 ft 

where xq = 15/xm and M = I0~ 26 kg (for lithium). While this is small compared to the 
phase shift which atom interferometric techniques p2|-p4| can obtain (~ 3 x 10 6 rad per 



g), improvements can be made. The size of the phase shift may be increased using more 
atoms, separating the wells further or by allowing a longer time r for the interaction. The 
number of atoms could be increased by up to 10 times without invalidating the two mode 
approximation. The two wells need to remain close during the preparation of the initial 
state and during tunneling, but could be separated and brought together againg during the 
weak force interaction. The time of the interaction is limited by mechanical vibration. 

Hence using a double condensate in a scheme to make sensitive measurements, such as 
that presented here, may be feasible. The experimental techniques currently being developed 
to produce and manipulate BECs may allow such a scheme to be realized in the near future. 
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FIGURES 



FIG. 1. Schematic outline of proposed measurement scheme. The stages are: (a) preparation 
of the superposition state, (b) interaction with the weak force for time r, (c) a ir/2 rotation caused 
by tunneling, and (d) a homodyne measurement of the number of atoms in one well. 

FIG. 2. The probability distribution resulting from the homodyne measurement. The fringes 
corresponding to adjacent m values can be distinguished if the coherent amplitude of the light field 
ao is large enough. 
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- homodyne measurement 




yCO^inxt/m-l) 



y(0)sin %tpn. y(0)sin Xtj(m+1) 



